Magnetization in two-dimensional electron gas in a perpendicular magnetic field: the 

roles of edge states and spin-orbit coupling 
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We study the de Haas-van Alphen (dHvA) oscillations in the magnetization of a two-dimensional 
electron gas (2DEG) under the influence of the edge states and/or the Rashba spin-orbit interaction 
(SOI). The boundaries of the systems lift partially the degeneracies of Landau levels (LL's) and the 
resulting edge states lead to the changes of both the center and the amplitude of the sawtoothlike 
magnetization oscillation. The SOI mixes the spin-up and spin-down states of neighboring LL's 
into two unequally spaced energy branches. The inclusion of SOI changes the well-defined sawtooth 
pattern of the dHvA oscillations in the magnetization. The weaker the magnetic field is, the larger 
is the change of the dHvA oscillations due to the edge effect and/or the spin-orbit coupling. Some 
theoretical results are compared with the experimental data. 

PACS numbers: 71.10.Ca, 73.20.At, 72.15.Gd 



CD 



I 

o 



> 
oo 

O 

o 



X 



INTRODUCTION 



The physics of two-dimensional electron gas (2DEG) 
in the presence of a perpendicular magnetic field reveals 
a rich variety of remarkable phenomena, for example the 
integer and fractional quantum Hall effects (QHEs). Re- 
cently, the magnetization property of 2DEG at low tem- 
perature and in a strong perpendicular magnetic field 
has attracted extensive interest. It is due to the fact 
that the magnetization is particularly suited to investi- 
gate the electronic ground-state properties and the den- 
sity of states (DOS) of 2DEG by its minimal pertur- 
bation to the system^. On the experimental side, the 
measurement of the magnetization keeps a most chal- 
lenging task due to weak signal of the magnetization. A 
variety of techniques, such as dc superconducting quan- 
tum interference device (SQUID) magnetometers^i^i^i^, 
picking up coils lithographed above the gateS, torque 



magnetometeri 



.7,8,9,10,11.12 



torsional magnetometer with 



optical angular detectiorti^iiiii^, and micromechanical 
cantilever magnetometer o"'^^i^^'^^i^^'^°i^^'^^ , have been de- 
veloped to observe the magnetization oscillations, i.e., 
the de Haas-van Alphen (dHvA) effect in high-mobility 
2DEG. Clear sawtooth dHvA oscillations in the mag- 
netization have been observed for LLs of filling fac- 
tors up to 53i^. More recently, a novel method has 
been used by Prus et al;^ and Shashkin et al.^^ to 
measure the spin magnetization of 2DEG in silicon 
metal-oxide-semiconductor field-effect transistors (MOS- 
FETs). This method entails modulating the magnetic 
field with an auxiliary coil and measuring the imagi- 
nary (out-of-phase) component of the ac current induced 
between the gate and the 2DEG system, which is pro- 
portional to d^/dB (where fi is the chemical potential 
and B is the magnetic field). Using the Maxwell rela- 
tion, djij dB~—dM. / dM , one can then obtain the mag- 
netization M. by integrating the induced current over 
the electron density, M . Pauli spin susceptibility has 



been observed to behave critically near the 2D metal- 
insulator transition, in agreement with previous trans- 
port measurements^i2^. With the similar method, Anis- 
simova et ali^ have studied the thermodynamic magne- 
tization of a low-disordered, strongly correlated 2DEG in 
silicon MOSFETs in perpendicular and tilted magnetic 
fields. By measuring dfi/dB at noninteger filling factors, 
they have directly determined the spectrum characteris- 
tics without any fitting procedures or parameters. 

On the theoretical side, extensive studies of dHvA os- 
cillations in the magnetization of 2DEG have also been 
carried ou t^^i^^i'^^i^^ '-^^''^-^--^^. In particular, Bremme et 
aL- - have investigated the influence of the edge current 
on the dHvA oscillations in the magnetization of a 2DEG 
using a spinless single-particle approach. Sharapov et 
al4i have extensively discussed the dHvA oscillations of 
the magnetization in planar system with the Dirac-like 
spectrum of quasiparticle excitation. In addition, the 
magnetization oscillations as a function of the magnetic 
field has also been theoretically studied in quantum dot 
systems22i^i2i. However, to our knowledge there are no 
detailed treatments of the influence of edge states and the 
SOI on the magnetization in 2D systems. In this paper, 
we study systematically the thermodynamic magnetiza- 
tion of a 2DEG system with edge states and SOI. In par- 
ticular, we address the effects of SOI and edge states on 
the Landau level (LL) structure, the chemical potential, 
and the magnetization and its susceptibility to strong 
magnetic field. Quantum oscillations in the magnetiza- 
tion of a 2DEG are well known to be characterized by 
strictly (l/i?)-periodic sawtoothlike oscillations with an 
amplitude of 1 effective Bohr magneton /is {—eh/ 2m* 
with m* the effective electron mass) per electron. We will 
show that the picture changes for the case in the presence 
of SOI and edge states. The degeneracy of LL's plays 
an important role in the formation of dHvA oscillations. 
The edge states lift partially the degeneracies of Landau 
levels and lead to the change of both the center and the 
amplitude of the sawtoothlike magnetization oscillation. 
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The SOI mixes the spin-up and spin-down states of neigh- 
boring LL's into two unequally spaced energy branches. 
The inclusion of SOI changes the well-defined sawtooth 
behavior of the dHvA oscillations in the magnetization. 
These results may be found useful in the characterization 
of magnetic oscillations in two dimensional systems. 

In Sec. II we review the exactly solvable cases of bulk 
2DEG with or without Rashba SOI, and numerically 
solvable cases in the presence of both edge states and 
Rashba SOI. In Sec. Ill we present the results for the 
magnetization and the effects of SOI and/or edge states. 
The results on the magnetic susceptibility are presented 
in section IV. Some concluding remarks are given in Sec. 
V. 



and the corresponding eigenstates are 



\n,k) = 
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), (3) 



I'^k is 



where Hn{x) is the Hermite polynomial, and 
the center of the cyclotron orbit. The degeneracy Nl of 
each Landau level per spin is given by Nl=2itS/11 with 
S=LxLy being the area of the sample. When the electron 
spin degrees of freedom is included but the edge-state 
effect and spin-orbit coupling are excluded, the eigenstate 
is characterized by three quantum numbers, 



|n, fc,s) = \n, k)xs 



(4) 



II. ENERGY SPECTRUM FOR 2DEG 

We consider a 2DEG with the Rashba coupling in the 
x-y plane of an area Lx x Ly subject to a perpendic- 
ular magnetic filed ^=Bz. The electrons are confined 
between and Ly in the y direction by an infinite po- 
tential wall, and its wave function is periodic along the 
X direction. We choose the Landau gauge A=—Byx. 
The Hamiltonian for a single electron of spin-1/2 with a 
Rashba coupling is given by 

^" = ■^—; + -^0^x<^y-^y'^x)^2^stJ^BB(T^ + V{y), (1) 

where m*, (— e), and gs are the electron's effective mass, 
charge and effective magnetic factor respectively, fiB is 
the Bohr magneton, Il=p + eA/c is the kinetic operator, 
A is the Rashba coupling, and CTq. are the Pauli matri- 
ces. The last term V{y) is the lateral confining potential: 
V{y)=0 tor ^ y ^ Ly and infinite otherwise. Rele- 
vant quantities related to the magnetic field are the cy- 
clotron frequency ujc—eB/m* and the magnetic length 
Ih^y/h/eB. The Rashba SOI in Eq. (1) stems fr om the 
structural inversion asymmetry (SIA) introduced by a 
heterojunction or by surface or external fields. In semi- 
conductors with narrower energy gap (InGaAs, AlGaAs), 
this effect is expected to be stronger—. It has been shown 
experimentally that the Rashba SOI can be modified up 
to 50% by external gate voltages^^'^''. 

Without considering the edge-state effect, the magne- 
tization of a 2D spinless electron moving in a high per- 
pendicular magnetic field is a 3D simplification that has 
been solved by Landau in a pioneering paper published in 
1930—. In that (2D) case, each state is described by two 
quantum numbers k and n. The quantum number k de- 
notes the X component of the electron momentum (scaled 
by h) and is a constant of motion. The other quantum 
number n being the Landau-level index describes differ- 
ent modes of a displaced linear oscillator of frequency ujc- 
The LL's are 



En 



{n + ^)hajc 



(2) 



where s = ±1/2 and x±i/2 is the eigenstate of spin oper- 
ator Sz with eigenvalues ±h/2. In this case, the spin-split 
LL's are 



(5) 



Thus each LL is split into spin-up and spin-down 
branches. When the edge channels are furthermore in- 
cluded but the SOI is exclude d ^"'" " , the eigenstate is 
still a product state of spin and orbital degrees of free- 
dom. However, the orbital part of the wave function is no 
longer a form of Hermite function due to the confinement 
along y direction. In this case, the eigenstate is given by 
|n,fc,s)=^e*''=^(p„,j^o(y)Xs, where (pn,yo{y) obeys 



2m* 



+ -m*u)c{y - yaf + V{y) 



(6) 



This equation has been solved by MacDonald et al. 
through properly applying the boundary condition**^. 
The resulting eigenvalue spectrum of Eq. ([6|) has a form 
^n{yo)=[t^n{yo) + \]^c, where Vniyo) is numerically ob- 
tained by requiring the wave function ^Pn,yo (j/) to vanish 
at the boundary ?;=0, Ly. In this case, the spin-split LL's 
are given by 



(7) 



which as an example is illustrated in Fig. 1(a) as a func- 
tion of the guiding center yo- 

On the other side, when the Rashba SOI is included 
while the edge-state effect is excluded, the spin-orbit cou- 
pling mixes the two spin components. In this case, the 
energies of the two branches of states (denoted by ±) are 
given by 



huj. 



(8) 



for n > I, where we have defined g=gsm* /2me and the 
effective (dimensionless) Rashba coupling 'q=\m* 1^ / . 
For n = 0, there is only one single state with the energy 
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EQ=fkJc{^ — g), which is the same as the lowest Lan- 
dau level without SOI. The corresponding eigenstates are 
given by 



l"'2^O'±>-(,^sin0±|n-l,yo; 



(9) 



for n > 1, where the parameters 9^ are given by 
tan6'^=— u„ ± y^l + with u„=(l — g)/\/^nri. For 
n=0, the single state is a product of the ground-state 
oscillator mode |0, fc) and eigenstate Xi/2 of Sz- Thus 
it is interesting to see that the ground state (n=0) has 
the fully-polarized spin along the z direction. In the ex- 
cited states the spin is tilted with an expectation value 



of its z component (ctz) =cos 



2 fl± . 



sin^ 9^ that decreases 



as A and n increase. A prominent feature is that the two 
branches of Landau levels and E~j^^ cross each other 
at the values of r] satisfying 

v/(l - 9? + 8n?72 + ^(1 - gY + 8(n + \)if ^ 2. (10) 

This degenerate behavior in the energy spectrum has 
been used to produce the resonant spin-Hall curren t^^'^^ . 

When both the edge states and the Rashba SOI are 
included, two interplayed mixing mechanisms occur. One 
is from the coupling between different confining orbital 
modes along the y direction; the other is the mixing of 
the eigenstates Xi/2 and X-i/2 of the spin operator s^. 
As a result, the wave function for the final Hamiltonian 
(1) can be written in a general form^^'^: 

1 



*n(a;,y) 



(11) 



Here fniu) are expanded in the basis of the infinite po- 
tential well, 




(12) 



with n being the Landau-level index and m an integer. 
The Schrodinger equation iJo^n=^^n(yo)^n leads to the 
following equations for the spinors: 



[Ai± - E„] 



air. 



(13) 



(14) 



Mirn = ^ (^^^ ^ d9sm{l9){9 - 0o)'sin(m6 



Elm. = 



'^^^^ f d9sin{l9){9-9a)sm{m9), 



7r2 I 



b Jo 



Gi^ ^ 27]^ f d9sin{l9)^sin{m9), 
Ly Jo 89 




\ ^mn J 

where cr* = {a^ ± i'^y)/'^: — En/fujJc- The other 
parameters in Eq. (13) are defined as 
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FIG. 1: The energy spectrum in units of meV versus the 
guiding center yo without (a) and with the Rashba spin-orbit 
couplings (b), respectively. In both figures, the parameters 
are chosen as m* = O.OSme, L = 600 nm, B = 1 T, and 
Qs — 4. The Rashba SOI strength in (b) is set as A = 15meV 
nm. 



where 9o—yaT^/Ly. We solve these equations in a trun- 
cated Hilbert space disregarding the states with energies 
higher than the cutoff energy. Typically we take a matrix 
Hamiltonian of dimension of a few hundred. We increase 
the size of the Hilbert space by a factor 2 and find no 
change in the results presented below. In all cases the 
width of the sample Ly is taken large enough to have 
the cyclotron radius Tc smaller than Ly/2. The right 
and left edge states are then well separated in real space. 
A typical energy spectrum is shown in Fig. 1(b). For 
yo — Ly/2 the states are equal to the bulk states, except 
for exponential corrections. The wave functions and the 
energy spectrum reproduces the above-discussed bulk re- 
sults without edge states. As yo approaches the sample 
edge, the effect of the confining potential becomes impor- 
tant and it generates the fc-dependcnt dispersion of the 
energy levels^, which has profound effects on magneto- 
transport and magnetization properties. 



III. THE MAGNETIZATION OF 2DEG 

We begin with a briefly review of the standard 
derivation of the 2DEG magnetization from the free 
energy. The magnetization density is the derivative 
of the Helmholtz free energy density with respect to 
B at fixed electron density Af and temperature T, 
M =~{dE/dB)\jyj-^T- For the present model with the LL's 
En.sivo), the Helmholtz free energy density is given by 

F{B, T) = ^AA - -L^ dyo J2^n{l + e/'[^-^"^-'(^o)] | 
P Jo t1 J 



M^-i J dED{E,B)\n{l + e^^''-^^, 



(15) 



where f3=l/kBT, N^=N]:^/S, and fj, is the chemical po- 
tential. Note that we have defined in the above equation 



the density of states (DOS) per area 

i^(£;,i3) = --i-^<5(i?-i?„,,(yo)) (16) 



LxLy — — 

^ n.k.s 



N, 



dyoS{E - En,s{yo)), 



where we have replaced fc-sum with J dk and used 
the relation yo=llk. The explicit inclusion of the DOS 
in the expression can be utilized to take into account the 
impurity effect, which broadens the LL's into Gaussian or 
Lorentzian in shape. For simplicity we did not consider 
the broadening effect in this paper. In the absence of 
edge states, the LL's -En.s(yo) Stve uniform in space and 
thus Eq. (fT5|) reduces to 



F(S,T) = ^AA-^^ln{l 



The B dependent chemical potential fi is connected to 
the experimentally accessible electron density Af via the 
local density of states (DOS). In the clean sample limit 
this is written as 



dyo ^ fns iVo) 



(18) 



where fns (vo)=^7rz — r-i — \ — 
for the spin-split LL's En.sivo)- 
netization density becomes 



is the Fermi distribution 
From Eq. (fT7)l the mag- 



n,s [ -^0 

+ r"^ln|l + e^[--^"^^-(-)l|" 

hfjjQ Ly I i 



(19) 



One can see that the magnetization consists of two parts. 
The first part M^^^^ is the conventional contribution from 
the B dependence of the LL's and thus denotes a diamag- 
netic response. The second part M^^^ comes from the 
B dependence of the level degeneracy factor N,y, thus 
describing the effect of the variation of the density of 
states upon the magnetic field and denoting a paramag- 
netic contribution to the total magnetization. Obviously, 
M(°) is negative while M^^^ is positive, the net result is 
an oscillation of the total magnetization M between the 
negative and positive values as a function of B. At zero 
temperature, the expression for M reduces to a sum over 
all occupied LL's: 



occ f „i 



dB 



(20) 
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FIG. 2: (Color online) Magnetic field dependence of (a) 
Chemical potential /x (thick line) and Landau levels (dashed 
line: spin down; solid line: spin up) and (b) magnetiza- 
tion m (per electron) for the 2DEG without edge states 
and spin-orbit coupling. The other system parameters are 
A/" = 4.5 X 10~''/nm^, = 4, m = O.OSrrie, for the inversion 
heterostructure Ino.53Gao.47As/Ino.52 Alo.48As, and r=3.5 K. 
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FIG. 3: (Color online) The chemical potential /i (a) and the 
magnetization per electron m (b) in units of /is in the sample 
with edges as a function of the external magnetic field B. The 
Rashba coupling A = 0. The sample size L — 600nm. The 
other parameters are the same as those in Fig. 2. 



where the sum runs over all occupied states and fiQ is the 
zero-temperature chemical potential (Fermi energy). 

For comparison, let us start from the conventional re- 
sult for the bulk 2DEG without SOI and edge-state ef- 
fects. In this case, both the chemical potential n [Fig. 
2(a)] and the magnetization (per electron) m [Fig. 2(b)] 
display the well-known sawtooth behavior with varying 
the magnetic field. At zero temperature, the explanation 
of the dHvA oscillation can be given with the help of 
the filling factor i'=Af/N^, which measures the number 
of occupied LL's and is an integer when all the available 
states in the v lowest Landau levels are filled. At these 
integer values the 2DEG is incompressible and the chemi- 
cal potential jumps discontinuously between two adjacent 
LL's by an amount of Afj,o—hLUc—2fiBB {fj,B—eB /2m*), 
which defines the incompressibility gap. Note that the 
abrupt jump in the dHvA oscillation is on the high mag- 
netic field side of the sawtooth, which is special for our 
present choice of the thermodynamic system. If the sys- 
tem is constrained to have constant chemical potential. 
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FIG. 4: (Color online) The magnetization per electron m in 
units of hb in the sample with edges as a function of the 
external magnetic field B. The Rashba coupling A = 0. The 
sample size L — 400, 600, 800nm. The other parameters are 
the same as those in Fig. 2. 

then the jump in the dHvA oscillation will move to the 
low magnetic field side of the sawtooth, which has been 
confirmed by Meinel et al.^ in an experiment with the 
electron density Af modulated by applying a gate volt- 
age to the sample. The zero-temperature behavior of the 
magnetization curve can be seen by Eq. (|20p . which in 
the absence of edge states is 

occu. 

M= ^ ^(Mo-2i;„^,). (21) 

n,s 

From above equation, one can easily derive a simplified 
Maxwell relation ■^^=^- Thus the discontinuous jump 
AM is related to the discontinuity in the chemical poten- 
tial A/i, i.e., the Landau energy gap at even and the Zee- 
man gap at odd filling factors, according to A/i = AM-B. 
When the magnetic field increases within an odd filling 
factor v, the Fermi energy and LL's also increase with 
B. As a result, the magnetization also increases rapidly 
as a function of B, then evolves a maximum at adja- 
cent filling factor ly — I and suddenly jumps to a nega- 
tive value. The zero-temperature magnetization jump at 
these even filling factors are given by the above Maxwell 
relation, AM=2/iB, independent of the magnetic field 
and spin-splitting. At finite temperature as shown in 
Fig. 2 (r=3.5 K), the oscillation amplitude of the mag- 
netization increases with increasing the magnetic field. 
This fact is due to that the influence of finite tempera- 
ture in this case (i.e., no SOI and no edge states in the 
2DEG) is merely to reduce the oscillation amplitude and 
the discontinuities in /i and M via the smearing of the 
Fermi-Dirac distribution. Another fact revealed in Fig. 
2 is that the inclusion of the Zeeman splitting in the LL's 
does not change the dHvA oscillation modes of the phys- 
ical quantities. It is due to this fact that the LL's are 
usually treated to be spin degenerate in previous work. 

Now let us see the edge-state effects on the magne- 
tization. Fig. 3 shows the influence of the edge states 
on the oscillations of chemical potential and magnetiza- 
tion (dHvA oscillations) with magnetic field. The most 
prominent feature brought by the edge states is that the 
center of the dHvA oscillations is now dependent on the 
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FIG. 5: (Color online) (a) The dependence of the center of 
dHvA oscillation on the size of the sample. The magntic 
field B is chosen aroud 1.2T. (b) The B dependence of the 
center of dHvA oscillations. In both figures the Rashba SOI 
is neglected. 




FIG. 6: (Color online) (a) Red curve: The bulk contribution 
of the magnetization m per electron (in units of /xs); Black 
curve: The total magnetization per electron m (in units of 
11b)- The Rashba coupling A = 0. (b) The ratio of edge states 
contribution of magnetization to the total magnetization. L = 
600nm. 



magnetic field. In particular, for the field less than 1 
Tesla, the oscillatory magnetization is always positive in 
sign. Another feature shown in Fig. 4 is that the os- 
cillation amplitude decreases with decreasing the sample 
size. As is known, the origin of dHvA oscillations is the 
degeneracy of Landau levels. The edge states with dis- 
persion lead to edge current which not only is crucial for 
the quantum Hall effects, but also very important for the 
magnetization'^'^. The dispersion of the edge states par- 
tially lift the degeneracy of the Landau levels. Thus the 
edge states tend to destroy the dHvA oscillations. There- 
fore it leads to the decreasing of oscillation amplitude as 
shown in Fig. 4. The upshift of the center of dHvA os- 
cillations may be understood in the following way. For 
the effects from the edge states, what really matters is 
the ratio of two important length scales: the magnetic 
length 4 and the system size Ly. The decreasing of Ly 
is equivalent to the increasing of Ij,, i.e., decreasing of 
B or ujc- As seen in Eq. (6), in the case with smaller 
Ly or weaker B, the second term becomes less impor- 
tant and the eigenenergy is less sensitive to the magnetic 
field for states with yo near the edges. From Eq. PU)) . 
one sees that the second term overcomes the first term 
and leads to the upshift of the center of dHvA oscilla- 
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tions. The smaller the system size is, the more profound 
effects the edge states lead to, as seen in Fig. 4 for both 
the center and amplitude of the dHvA oscillations. Fig. 
5(a) shows quantitatively the system size dependence of 
the shift of the oscillation center. It has the dependence 
Roughly, the contribution of the edge states is 
proportional to the number of edge states (as also seen 
from Eq. (|20p ). which is proportional to vTc/Ly, where 
the cyclotron radius Tc = -Jvlb^, with the number of 
the occupied Landau levels v ~ Thus the center 

of dHvA oscillations is proportional to l^/Ly — 1/B^Ly. 
The B and Ly dependence is clearly seen in Figs. 4 and 5. 
To see more explicitly the contribution from edge states 
and bulk states, we plot the total magnetization and the 
contribution from bulk states in Fig. 6(a). The contri- 
bution from the edge states is obtained from Eq. (fT5)) by 
summing over terms from edge states, with |?/o| < fc or 
\L — yo\ < Tc- The rest contribution is from bulk states. 
There is no upshift of the magnetization oscillation cen- 
ter for the part from bulk states. It shows explicitly that 
the upshift of the center of dHvA oscillations is due to 
the existence of edge states. Fig. 6(b) shows the depen- 
dence of edge states contribution on the magnetic field. 
The contribution from edge states increases as decreasing 
the magnetic field, or equivalently decreasing the sample 
size as one expects. 

When the Rashba SOI is introduced, the filling factor i' 
is not linearly proportional to the inverse of the external 
field B, and there is a energy competition between the 
Zeeman coupling and the Rashba coupling. Also, due to 
the entanglement between the orbital and spin degrees of 
freedom, it is difficult to distinguish their separate contri- 
butions to the total magnetization. These factors make 
the physical picture of the dHvA oscillations to change 
fundamentally, as shown in Fig. 7(a) for chemical poten- 
tial fi and Fig. 7(b) for magnetization M as functions of 
B. One can see from these two figures that the Rashba 
SOI has no visible influence on the magnetic oscillations 
of the quantities /i and M at large values of B, where 
the Zeeman and spin-orbit coupling splitting are small 
compared to the Landau level splitting. At low magnetic 
field, however, the SOI modulation of the magnetic os- 
cillations becomes obvious, which can be clearly seen by 
the enlarged plots of /i and M in Fig. 7(c) and 7(d) 
respectively for B less than 2.4 T. For comparison, we 
also re-plot in Figs. 7(c)-(d) the cases without Rashba 
SOI. One can see from these two figures that the SOI 
brings about two new features at low magnetic field: (i) 
The sawtooth structure is inversed, i.e., the location of 
peaks in ^ and M with SOI correspond to the valleys 
without SOI. This inversion is due to the different LL's 
in the two cases, (ii) The oscillation mode is prominently 
modulated by SOI and a beating pattern appears. This 
beating behavior in the oscillations are due to the fact 
that the LL's and E~ are now unequally spaced due 
to the presence of SOI. 

Another effect caused by the Rashba spin-orbit cou- 
pling is that there are weak peaks appearing in the chem- 
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FIG. 7: (Color online) (a) The Chemical potential fi and (b) 
the magnetization m (per electron) for the 2DEG under dif- 
ferent temperatures T — 3.5K and IK with Rasha spin-orbit 
interaction in the absence of edge states, the Rashba cou- 
pling is taken to be A=15 meV nm. The other parameters 
are same as those in Fig. 2. The detailed oscillations of fi 
and m at low magnetic fields are illustrated in (c) and (d) 
respectively (solid lines), in comparison with the case of A=0 
(dotted fines). 
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FIG. 8; (Color online) The chemical potential fi (a,c) and 
the magnetization m (per electron in units of hb) (b,d) in 
the sample with edges as a function of the external magnetic 
field B. The Rashba coupling A = (a,b) and A = 15meV 
nm (c,d). In four figures, the sample size L = 600nm. The 
temperature T — 3.5 K. 



ical potential and the magnetization oscillation at low 
magnetic field, as shown with the arrows in Fig. 7. In 
the 2DEG sample without edges, the weak peaks will 
appear in the magnetization once the temperature is suf- 
ficient low, for example, T=1K. These weak peaks have 
been observed recently by Schaapman et alJ^ when they 
measured the magnetization of a dual-subband 2DEG, 
confined in a GaAs/AlGaAs heterojection, and by Zhu 
et al.^*^ when they measured the magnetization of high- 
mobility 2DEG. These peaks is so weak that they will 
disappear when the temperature is much higher. 

Now let us see the SOI effect superposed on the edge- 
state effect. Figures 8(a)-(b) show the chemical potential 
fj, and magnetization M respectively in the presence of 
both SOI and edge states. One can see from these two 
figures that at large values of B, where the Zeeman and 
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spin-orbit coupling splitting are small compared to the 
Landau level splitting, and thus the Rashba SOI has no 
big influence on the magnetic oscillations of the quanti- 
ties /i and M . Neither do the edge states, since the cy- 
clotron radius is much smaller than the system size. At 
low magnetic field, however, both edge states and Rashba 
SOI change the pattern of dHvA oscillations, as clearly 
seen in Figs. 8(a), 8(b), and the enlarged plots of /j, and 
M in Figs. 8(c) and (d) respectively for B less than 2.4 
T. For comparison, we also re-plot in Figs. 8(c)-(d) the 
cases without Rashba SOI. One can see that the total ef- 
fects are the superposition of effects from both the edge 
states and Rashba SOI. 



IV. MAGNETIC SUSCEPTIBILITY OF 2DEG 

Now we turn to study the magnetic susceptibility x{B) 
for the 2DEG. From Eq. (fT9|l . one obtains the expression 
for x(^) as follows: 



(22) 




f3 



2 /3[i;„,efao)-A'] 



' dB 
dB 



+/«s(yo) 



which at zero temperature reduces to 



xiB) 



occu 




^ dB 



(23) 



For the ideal noninteracting 2DEG without the edge 
states, the second term in Eq. is zero due to the 

fact that both the Landau and the Zeeman splittings of 
the energy spectrum are linear in B. Thus in this case 
(zero temperature and no edge states) the magnetic sus- 
ceptibility is simply written as 



xiB) 



occu 

e X ^ 



h ^ [dB 



,dE^ 
' dB 



(24) 



Furthermore, if the SOI is disregarded, the derivative 
Oiiq/OB in quantizing magnetic fields (except at even 
integer filling factors) is equal to 



dfiQ 
dB 



= Ms 



1 , 2mp 
2 m* 



1 

T^9s 



(25) 



Then Eq. (IM)) clearly shows Landau diamagnetic (x < 0) 
at Landau gaps and Pauli paramagnetic (x > 0) re- 
sponses at Zeeman gaps. 



J3" 
CD 




J2 




, (a) 

0.2 OA 0.6 0.8 1.0 
1/B(1/T) 



(b) 

0.4 0.6 0.8 1.0 
1/B(1/T) 



FIG. 9: (Color online) The magnetic susceptibility x in the 
sample (a) without and (b) mth edges as a function of the 
inversed magnetic field 1/B. The black (real) line is cor- 
responding to the case of the Rashba SOI strength A = 0, 
while the red (dotted) line is corresponding to the case of 
A = 15meV nm. In both figures, the electron density and the 
temperature are respectively set as Us = 4.5 x 10~^/nm^ and 
T = 3.5K. The system size in (b) L = 600nm. 



Let us first see the SOI effect on the magnetic suscep- 
tibility. Figure 9(a) shows x [B] as a function of 1/ B for 
the edgeless 2DEG. One can see from Fig. 9(a) that there 
are a series of equal-distance resonance peaks appearing 
in the magnetic susceptibility x with the magnetic field 
with/without the Rashba coupling A. With the magnetic 
filed increasing, the magnitude of the resonance suscep- 
tibility increases. Similar to that of the magnetization, 
the explanation of the resonance peaks of the magnetic 
susceptibility without Rashba coupling also needs the 
help of the filling factor v ^ 1/B. When v increases 
to an integer value, the 2DEG is incompressible. At this 
time, all the available states in the v lowest Landau lev- 
els are filled. Upon increasing the inverse magnetic field 
1/B, the electrons are transferred to the next Landau 
level. Thus the chemical potential changes discontinu- 
ously (See Fig. 2(a)) and the magnetization jumps (See 
Fig. 2(b)). Therefore there is a corresponding resonance 
peak appearing in the magnetic susceptibility x- When 
the Rashba spin-orbit coupling is introduced, it will de- 
stroy the simply linear relation between the filling factor 
V and the inverse magnetic field 1/B. When the value of 
1/B is larger, the SOI effect is more evident (See the red 
line in the Fig. 9(a)). 

Figure 9(b) plots the magnetic susceptibility x in a 
2DEG sample with edges as a function of the inversed 
magnetic field 1/B. Comparing the black lines in Figs. 
9(a) and 9(b), which corresponding to the case without 
Rashba spin-orbit coupling, one can find that the reso- 
nant (for magnetic susceptibility) magnetic field shifts to 
larger values for the case with edges. We have obtained 
the conclusion that the edges not only affect the center 
of the dHvA oscillation, but also the values of the mag- 
netic field when the magnetization has the discontinuous 
change. As a consequence, the conditions of the suscep- 
tibility (defined as x=dM / dB) having resonance are the 
same with those of the magnetization. When the Rashba 
spin-orbit coupling is introduced, similar to the case in 
the sample without edges, the Rashba coupling changes 
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the resonance condition. When the vahie of 1 /B is much 
larger, the change is much larger as seen in Fig. 9(b). 

V. SUMMARY 

In summary, we have systematically studied the dHvA 
oscillations of the magnetization and its susceptibility for 
the 2DEG with the edges states and SOI included in the 
system. 

We find that the edge states and Rashba SOI play im- 
portant roles when the external magnetic field is small. 
The edge effect prominently changes the oscillation cen- 



ter and oscillation amplitude. The dHvA oscillation will 
change the sawtoothlike form if the Rashba coupling in- 
troduced, no matter the sample is with or without edges. 
The total effects are the superposition of effects from 
both edge states and SOI. 

Acknowledgments 

This work was supported by NSFC under Grants 
No. 10604010, No. 60776063, No. 10744004, and No. 

10874020, and by the National Basic Research Program 
of China (973 Program) under Grant No. 2009CB929103. 



* Author to whom correspondence should be addressed. 
Email address; zhang_ping@iapcm.ac.cn 

^ D. Shoenberg, Magnetic Oscillations in Metals (Cam- 
bridge University Press, Cambridge, 1984). 

^ H. Stormer, T. Haavasoja, V. Narayanamurti, A. C. Gos- 
sard, and W. Wiegmann, J. Vac. Sci. Technol. B 1, 423 
(1983). 

^ I. Meinel, D. Grundler, S. Bargstadt-Franke, Ch. Heyn, 
and D. Heitmann, Appl. Phys. Lett. 70, 3305 (1997). 

"* 1. Meinel, T. Hengstmann, D. Grundler, D. Heitmann, 
W. Wegscheider and M. Bichler, Phys. Rev. Lett. 82, 819 
(1999). 

^ I. Meinel, D. Grundler, D. Heitmann, A. Manolescu, V. 

Gudmundsson, W. Wegscheider, and M. Bichler, Phys. 

Rev. B 64, 121306(R) (2001). 
® F. F. Fang and P. J. Stiles, Phys. Rev. B 28, 6992 (1983). 

J. P. Eisenstein, Appl. Phys. Lett. 46, 695 (1985). 
® J. P. Eisenstein, H. L. Stormer, V. Narayanamurti, A. Y. 

Cho, A. C. Gossard, and C. W. Tu, Phys. Rev. Lett. 55, 

875 (1985). 

^ S. A. J. Wiegers, M. Specht, L. P. Levy, M. Y. Simmons, 
D. A. Ritchie, A. Cavanna, B. Etienne, G. Martinez, and 
P. Wyder, Phys. Rev. Lett. 79, 3238 (1997). 

^° M. Zhu, A. Usher, A. ,J. Matthews, A. Potts, M. Elliott, W. 
G. Herrenden-Harker, D. A. Ritchie, and M. Y. Simmons, 
Piiys. Rev. B 67, 155329 (2003). 

" A. J. Matthews, K. V. Kavokin, A. Usher, M. E. Port- 
noi, M. Zhu, J. D. Gething, M. Elliott, W. G. Herrenden- 
Harker, K. PhilUps, D. A. Ritchie, M. Y. Simmons, C. 
B. Sorensen, O. P. Hansen, O. A. Mironov, M. Myronov, 
D. R. Leadley, and M. Henini, Phys. Rev. B 70, 075317 

(2004) . 

D. R. Faulhaber and H. W. Jiang, Phys. Rev. B 72, 233308 

(2005) . 

M. R. Schaapman, P. C. M. Christianen, J. C. Maan, 
D. Reuter, and A. D. Wieck, Appl. Phys. Lett. 81, 1041 

(2002). 

^'^ M. R. Schaapman, U. Zeitler, P. C. M. Christianen, and 
J. C. Maan, D. Reuter, A. D. Wieck, D. Schuh, and M. 
Bichler, Phys. Rev. B 68, 193308 (2003). 
I. M. A Bominaar-Silkens, M. R. Schaapman, U. Zeitler, 
P. C. M. Christianen, J. C. Maan, D. Reuter, A. D. Wieck, 
D. Schuh, and M. Bichler, New J. Phys. 8, 315 (2006). 
M. P. Schwarz, D. Grundler, I. Meinel, Ch. Heyn, and D. 
Heitmann, Appl. Phys. Lett. 76, 3564 (2000). 



" M. P. Schwarz, M. A. Wilde, S. Groth, D. Grundler, Ch. 
Heyn, and D. Heitmann, Phys. Rev. B 65, 245315 (2002). 

^® M. A. Wilde, M. P. Schwarz, Ch. Heyn, D. Heitmann, and 
D. Grundler, D. Reuter and A. D. Wieck, Phys. Rev. B 
73, 125325 (2006). 

1" R. Knobel, N. Samarth, J. G. E. Harris, and D. D. 
Awschalom, Phys. Rev. B 65, 235327 (2002). 

2" J. G. E. Harris, R. Knobel, K. D. Maranowski, A.C. Gos- 
sard, N. Samarth, and D. D. Awschalom, Phys. Rev. Lett. 
86, 4644 (2001). 

M. A. Wilde, M. Rhode, Ch. Heyn, D. Heitmann, D. 

Grundler, U. Zeitler, F. Schaffler, R. J. Haug, Phys. Rev. 

B 72, 165429 (2005). 
22 N. Ruhe, J. I. Springborn, Ch. Heyn, M. A. Wilde, and D. 

Grundler, Phys. Rev. B 74, 235326 (2006). 
2^ O. Prus, Y. Yaish, M. Reznikov, U. Sivan, and V. Pudalov, 

Phys. Rev. B 67, 205407 (2003). 
^'^ A. A. Shashkin, S. Ariissiniova, M. R. Sakr, S. V. 

Kravchenko, V. T. Dolgopolov, and T. M. Klapwijk, Phys. 

Rev. Lett. 96, 036403 (2006). 
2^ S. V. Kravchenko and M. P. Sarachik, Rep. Prog. Phys. 

67, 1 (2004). 
2« A. A. Shashkin, Phys. Usp. 48, 129 (2005). 

S. Anissimova, A. Venkatesan, A. A. Shashkin, M. R. Sakr, 

and S. V. Kravchenko, and T. M. Klapwijk, Phys. Rev. 

Lett. 96, 046409 (2006). 
2« D. Shoenberg, J. Low Temp. Phys. 56, 417 (1984). 
2^ F. Oppen, D. UUmo, and H.Baranger, Phys. Rev. B 62, 

1935 (2000). 

3° L. Bremme, T. Ihn, and K. Ensslin, Phys. Rev. B 59, 7305 

(1999). 

31 S. G. Sharapov, V. P. Gusynin, and H. Beck, Phys. Rev. 
B 69, 075104 (2004). 

32 A. Aldea, V. Moldoveanu, M. Nita, A. Manolescu, V. 
Gudmundsson, and B. Tanatar, Phys. Rev. B 67, 035324 
(2003). 

33 W.-C.Tan and J. C. Inkson, Phys. Rev. B 60, 5626 (1999). 
3'' U. Sivan and Y. Imry, Phys. Rev. Lett. 61, 1001 (1988). 
35 E. I. Rashba, Sov. Phys. Solid State 2, 1109 (1960). 

3« J. B. Miller et al, Phys. Rev. Lett. 90, 076807 (2003). 
3'' J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. 

Rev. Lett. 78, 1335 (1997). 
3* L. Landau, Z. Phys. 64, 629 (1930). 
3'^ B. I. Halperin, Phys. Rev. B 25, 2185 (1982). 

A. H. MacDonald and P. Streda, Phys. Rev. B 29, 1616 



9 



(1984). 

S.-Q. Shen, M. Ma, X. C. Xio, and F. C. Zhang, Phys. 
Rev. Lett. 92, 256603 (2004). 

Z. Wang and P. Zhang, Phys. Rev. B 75, 233306 (2007). 



A. Rcynoso, G. Usaj, M. ,J. Sanchez, and C. A. Balseiro, 
Phys. Rev. B 70, 235344 (2004). 



